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Theory of a weakly non-ideal Bose gas in the canonical ensemble is developed 
without assumption of the C-number representation of the creation and annihilation 
operators with zero momentum. It is shown that the pole of the "density-density" 



1-^ ' Green's function exactly coincides with the Bogolybov's phonon-roton spectrum of 

Q . excitations. At the same time, a gap exists in the one-particle excitation spectrum . 



This gap is related to the density of particles in the " condensate" . 
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INTRODUCTION 



^ ! Starting with the Bogolyubov's papers [1,2], the microscopic theory of the degenerate 

H ■ 
.^.' Bose gas has been based on the special assumption that the creation a^ and annihilation ag 

operators of particles with zero momentum can be replaced by a C-number 

a+ = ao = (< A^o >)^/2, (1) 

where < A'q > is the average quantity of particles in the state with p = ("condensate"). 

This assumption leads to the necessity of introducing in the theory the anomalous aver- 
ages ("quasi- aver ages"), which is unnatural for the homogeneous and isotropic system under 
consideration. Let us pay attention that statement (1) in Bogolyubov's papers is directly 



related to the requirement of the weak interparticle interaction in the degenerated Bose gas. 
Though the Bogolyubov theory (see, e.g., [3,4]) gives rise to many important and widely 
recognized results, such as the expression for the spectrum of excitations, explanation of 
the experimental data and agreement with the Landau superfluidity condition, there are 
serious doubts on the validity of the relation as well as on the agreement of the Hamilto- 
nian corresponding to the assumption ([1]) with the original Hamiltonian of the system under 
consideration [5,6]. 

In this work, we propose a selfconsistent noncontradictory description of the degenerate 
Bose gas in which we do not use the assumption ([T]). In this way we show, in particular, 
that there exists the gap in the spectrum of one-particle excitations. The statement that the 
spectrum of one-particle excitations has a gap in parallel with the usual phonon-roton branch 
of the collective excitations in the degenerate Bose gas has been formulated in [7] . Later, 
the existence of the gap has been suggested in [8] . Spectra, thermodynamics, and dynamical 
structure factor for weakly non-ideal Bose gas below the condensation temperature T <Tq 
have been considered in [9-11] in terms of the dielectric formalism, similar to that in plasma- 
like systems. The dielectric formalism based on the Bogolyubov assumption ([T]) for the 
operators a^ and a^ was developed in [12]. Recently, the gap-existence problem for the 
weakly non-ideal Bose gas has been considered by different methods in papers [13-16]. In 
fact, the paper [14] confirmed (independently) and developed not only the statement of 
Ref. [7] about the existence of a one-particle gap in the weakly non-ideal Bose gas below 
the condensation temperature, but also the results of the papers [9-11] on the possibility of 
considering the superfiuid system without symmetry breaking. Surely, the existence of the 
gap disagrees with the conventional opinion that the gap is missing [3]. But we do believe 
that the conclusion about a missing gap is related to the Bogolubov simplifying assumption 
dl]), and is not obliged to persist in a more rigorous theory. 

It is necessary to stress that in the previous papers both Landau and Bogolyubov admitted 
that the one-particle excitation spectrum can have a gap, but later they dropped this idea, 
because of contradiction with observations of the phonon-roton branch of spectrum at small 
q in the neutron scattering experiments. The approach which includes the both spectra, with 
and without a gap, has not been found and coexistence of these two branches of excitations 
was not suggested. 

In the present work, we suggest that in the limit of strong degeneracy T — t- 0, where T is 



the system temperature, all the particles tend to occupy the zero-momentum state, which 
means that 

lim < No >=< N >, (2) 

where < A^ > is the average number of particles in the system. This statement is confirmed 
in the paper by the self-consistent consideration. Below we use the canonical ensemble where 
A^ is a given C-number. 

II. STATISTICAL SUM AND AVERAGES IN CANONICAL ENSEMBLE 

Let us consider the Hamiltonian of the non-ideal Bose gas of particles with zero spin and 
a mass m in a volume V 

if = ^£pa+ap + — J2 «(9)«pi-q/2«p2+q/2«P2-q/2«pi+q/2, (3) 

p q,pi,P2 

where aj" and ap are the creation and annihilation operators of particles with the momentum 

hp, 

[ap2,ap^J = ctpjO-p^ — '3'p;^Q'p2 = 'Jpi,P25 [^) 

Ep = h p^/2m is the energy spectrum of a free particle and u{q) is the Fourier-component 
of the inter-particle interaction potential. 

It is convenient to extract the particular term Uq with q = from the sum over q in the 
Hamiltonian (El). This term can be written as 



Uo = m(0)^^^, N = J]a+ap, u{0) = u{q = 0) = u{q ^ 0), (5) 

p 

where A^ is the operator of the total number of particles, and N is the C-number N =< N >. 

Hereafter, the brackets < ... > mean the canonical-ensemble averaging. Since Uq is also the 

C-number, one can write the statistical sum of the system under consideration in the form 

Z = Sp exp{-H/T) = exp (-^0) ^^^^" ^^ ) Sp exp(-i7o/T), (6) 

where 

Ho = H-Uo = J2^p(^p(^P + ^ Yl ^(^)"Ji-q/2«j2+q/2«P2-q/2«pi+q/2' (7) 

p q7^0,pi,p2 



Therefore, to provide convergence of the statistical sum ([7]) in the thermodynamic hmit 
y — 7- oo A^ — 7- oo, n = N/V = const., the known condition m(0) > has to be fulfilled. 
As f/o is a C-number, it does not affect any averaging at all, and for calculating of average 
values the Hamiltonian H is equivalent to Hq. For an arbitrary operator < A >, we get 

<A>= Z-^Sp {exp{-H/T)A} = Z^^Sp {exp{-Ho/T)A} =< A >o; Zo = Sp exp{-Ho/Tl8) 

Let us now consider a more complicated situation associated with calculation of a time 
dependent correlation functions f{t) such as 

f{t) =< [A{t),B{0)] >, A{t) = exp{iHt/h)Aexp{-iHt/h). (9) 

The time dependence of the operators flp (t) and ap{t) can be represented in the form 

a+(t) = exp{iHot/h)a^ exp{-iHot/h) exp{iNu{0)t/Vh). (10) 

ap{t) = exp{-iNu{0)t/Vh) exp{iHot/h)apexp{-iHot/h). (11) 

Therefore, if in each of the operators A and B the numbers of creation and annihilation 
operators coincide (which is typical of the operators of the physical variables), the time- 
dependent correlation function can be written as 

fit) =< [A{t),B{0)] >o, A{t) = exp{iHot/n)Aexp{-iHot/n). (12) 

On this basis, in what follows, we consider the average values with the Hamiltonian Hq (C]) 
within the canonical ensemble. The free energy of the initial system with the Hamiltonian 
H ([3]), according to ©-(J?]), reads 

F = -TlnZ = Uo + Fo, Fq = -TIuZq. (13) 

III. EQUATIONS FOR "DENSITY-DENSITY" GREEN FUNCTION 

Experimentally, the spectrum of collective excitations is found usually from data on the 
well observable maxima in the dynamical structure factor ^(q, w) for q 7^ 0, 

1 /"^ 

Si(l,uj) = 77 / expiiut) < p^{t)p-q{0) >o dt, (14) 

^ J —00 



PqW = Yl ^p-q/2(^)«P+q/2(^). (15) 

P 

where Pq(t) is the Fourier- component of the density operator in the Heisenberg representa- 
tion. The dynamical structure factor ^(q, u) ( !T^ is directly related to [17] with the retarded 
density- density Green's function x^(Q; ^) which is analytical in the upper semi-plane of the 
complex variable z {Imz > 0), 

2h 

X^(q, Z) = -— / dt exp{izt) < [Pq(t)p_q(0)] >0= - « Pq | P-q »„ (17) 

The definitions (TT6|) . ( TT71) have to be taken in the thermodynamic limit, where 

lim < No >o=< N >o= A^. (18) 

According to Eq. flT4l) . the retarded function x^{li ^) can be represented in the form 



^' ^) = t7 5^1 ^^P' ^' ^)' ^(P' ^' ^) = ^^ «p-q/2«P+q/2 I P-q))z (19) 



P 



The equation of motion for the function -F(p, q, z) with the Hamiltonian Hq, determined by 
(Eq. ([7])), can be written in the form 

{hz + £p-q/2 - £p+q/2) i^(p, Q, 2:) = /p_q/2 " /p+q/2 " 

— 2_]m(A;) 2_](( (ap+k-q/2«pi-k/2"pi+k/2«P+q/2 " «p-q/2^pi-k/2^Pl+k/2«p-k+q/2) | P-q)l2'^) 
k pi 

Here /p is the one-particle distribution function over the momenta ^p 

/p =< a+ap >o (21) 

For a temperatures T < Tq, where Tq is the condensation temperature, the one-particle 
distribution function /p can be represented as [18] 

/p=<Aro>Vo + /J(l-Vo), (22) 

where A'o = a,|ao is the operator of the number of particles with zero momentum ("conden- 
sate") , fp =< ctpCtp >o is the one-particle distribution function with non-zero momenta 
(the "overcondensate" states). Therefore, 

- = -o + ^E/p=-o+/7C?^/J, (23) 






where no =< Nq > jV is the average density of particles in the condensate. From ( 120|) 
and fl22l) . we find that the function -F(p, q, z) has singularities at p = ±q/2. Therefore, the 
density-density function x^(Q'; A (fT9ll can be written as 

X\q. z) = ^ F{q/2, ^,z) +^ F(-q/2, q,z) +^ J^ ^"^(P' "^^ ^) (^4) 

P7^±g/2 

The index T means, that the respective function describes the " overcondensate" particles. 
The singularities (conditioned by the condensate) in this function are absent. Then in the 
last term in Eq. fl24|) . we can change summation by integration over momenta. The functions 
F(±q/2, q, 2), extracted above, satisfy, according to Eqs. (!20|) .( l22|) . to the exact equations 
of motion 



V 



{hz - Eq) F(q/2, q, z) = [{Nq) - /^] - ^J2^^^^ J2^^ (4«pi~k/2«pi+k/2aq - 

pi 

«C^«pl-k/2«Pl+k/2aq-k) I P-q))z (25) 

{hz + Sq) F(-q/2, q, z) = -[{No) " /q ] - ^7 5Z "(^^ ^^^ (4-q«p,-k/2«pi+k/2ao - 

«-q«Jl-k/2«Pl+k/2a-k) I P-q))z- (26) 



V 

k^O pi 



Let us further consider the case of strongly degenerate gas, where T — )■ 0. To find for this 
case the main terms, following to the Bogolyubov's procedure, let us extract the terms on 
the right-hand sides of Eqs. ( !25l) .( !26l) . which are determined by the maximum number of 
the operators a^ and oq. In the limit of strong degeneration, taking into account ( IT8l) . we 
can omit the other terms in the course of calculation of F. Then, since f^ = f^ (q 7^ 0), 
Eqs. ([25D,(I26D take the form 

{hz T Sq) F(°)(±q/2, q, z) = ±[{No) - /J] ± ^u{q){{ (a+a+a^ao + a+al^aoao) | p^q))i27) 

Eqs. f l27|) are exact at low temperature, since they define the functions F^^'{±q/2,q^, z). 
For calculation the Green's functions in the right side of Eqs. (p71) . it is necessary to make 
some approximations. 

IV. DETERMINATION OF THE COLLECTIVE EXCITATIONS 

It should be emphasized that the Bogolyubov's approach, in which the operators a^ and 
Qq are considered as C- numbers, leads to violation of the exact relations (l25l) and (|26|) . In 



this approach ([T]), the main term F^^^ of the function F has the form [19] 

F(°)(±q/2,q,^) = {No) ((a±q | a+q)), (28) 

Respectively, instead of the exact equations for two-particle Green's functions (l25l) and ( l26l) . 
we obtain the equation of motion for one-particle Green's functions form (l28l) . As is shown 
below, the equations of motion for the two-particle Green's functions without the approx- 
imation on the C-number representation of the operators a^ Oq are essentially different 
from the equations for the one-particle distribution functions. Therefore, for calculation of 
the Green's functions in the right part of ( !27|) . the assumption ([T]) cannot be used. At the 
same time, the idea of the C-number approximation for some operators seems, itself, very 
attractive. Below we apply this idea to the varsion alternative to the Bogolyubov assump- 
tion. According to flTSj) . it is natural to accept that to calculate the Green's functions in 
the limit of strong degeneracy T — )■ 0, the C-number approximation has to be applied not 
to the operators a^ ag, but to the physical value - the operator of the number of particles 
in the "condensate" A^o 

No = {No). (29) 

For low temperatures T, the operator A"o in the averages (under the angle brackets) can 
be changed by the total density operator A^ (equal to Pq=o), which is the C-number in 
the canonical ensemble. Therefore, in the case under consideration, it can be removed from 
brackets. Therefore, for low temperatures, the proposed approximation seems asymptotically 
exact. 

In this approximation from (127|) directly follows 

{hz - e,) FW(q/2, q, z) = [{No) - /J] + ^^n(g){FW(q/2, q, z) + F^'\-q/2, q, zJ^Q) 

{hz + e,) FW(-q/2, q, z) = ~[{No) - /J] - ^^«(g){FW(q/2, q, z) + FW(-q/2, q, zj^l) 

From Eqs. (130|) . (!3T]) . one can find the solutions for the functions F{±q/2, q, z) 

m 10 f,\ [<N,>-fl]{nz + e,) [<No>-fl]{hz-e,) 



huj{q) = ^el + 2nou{q)eq (33) 

The relation (|33ll for the spectrum hu{q) corresponds exactly to the known Bogolyubov 
expression [1,2]. By substituting to ( TT9l) . and taking into account that for the case of 
strong degeneration the contribution of the functions -F^(p, q, z) is negligible, we obtain the 
expression for the main term x^^Kq^ -2) of the "density- density" Green function x(g, z) 

As it is well known, the singularities of the function x^(Q') ^) determine the spectrum of 
collective excitations in the system. Therefore, under the assumption about C-number 
behavior of the operator A^q; w6 obtain the Bogolyubov's result for the spectrum of the 
collective excitations in the degenerated and weakly interacting Bose gas. However, the 
question on the braced term f^/ < Nq > in flM|) remains open. The problem is in the 
behavior of the function f^ for the ideal Bose gas [20] 

r:-U('-^)-iV (35) 



T 



In the limit of small wave vectors q the function /* converges at non-zero temperatures 



q 



(hereafter, 1/g -divergence). Moreover, 



lim lim P'^ ^ lim lim P*^ (36) 

The similar problem arises when ( 122|) is used. 

V. ONE-PARTICLE EXCITATIONS 

To calculate the distribution function fj for the " overcondensate" particles in Bose gas, 
let us consider the one-particle Green function g^{q, z) 

(7^(g,z) = ((aq|a+))„ q^O (37) 

This function is directly connected with the distribution function /J by the relation [21] 

fj = j" ^g''{q,u;), g<{q,u) = -2hLxp(^]-l\ Img''{q,u + zO) (38) 



The equation of motion for the Green function (7^(g, z) for q ^ reads 

{hz - Eg) g^{q, z) = l + -Y^ u{k) 5]](( flp+kapflq+k | a+)),. (39) 

kjtO p 

As for the " density-density" Green function, we consider the case of a strong degeneration 
and extract the terms with maximum quantity of the operators a^ and Oq in the right-hand 
side of Eq. (139|) . Then from Eq. (139|) , we find 

{hz-Sq) g^{q,z) = 1 + —u{q) {{{a^aq^ao \ aq))^ + {{at^aoao \ a+))J (1 -5q,o), (40) 

In this case, it is assumed that the separated main term on the right-hand side of (40) 
corresponds to the case of the weak interaction. In the general case, when the interparticle 
interaction potential u{q) is not small, in separating the terms with a maximum number of 
the operators a^ and a^ on the right-hand side of (40), it is necessary to write an infinite 
series on the perturbation theory on the interaction u{q). It is clear that a similar situation 
occurs when using (27) instead of (26). Thereby, strictly speaking, the question whether the 
terms separated on the right-hand side of (40) are main from the viewpoint of the passage 
to the limit of the weak interaction remains open. 

If we now apply the assumption (1) about the C- number representation of particle cre- 
ation aQ and annihilation oq operators at a zero momentum to the calculation of Green's 
functions on the right-hand side of (40), we come to the necessity of introducing the so-called 
"anomalous" Green's functions 

9anomil, ^) = (( ^Iq | <))., Q ^ 0, (41) 

Let us pay attention that the nonzero value of the anomalous function Qanomi strictly speak- 
ing, requires the influence of an external field of special form on the system under consider- 
ation. As a result, taking into account (1) and (41), relation (40) takes the form 

{hz - Eq) g^{q, z) = l + nou{q) {g^{q, z) + g^^,^{q, z) ] . (42) 

Then we use the equation of motion for the anomalous function ganomilj z) at q y^ with 
Hamiltonian (7), 

{hz + Eq) ganomil^z) = -y^^u{k){{a^_^^a:^_j^ap I a+))„ (43) 
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Now, as in the above consideration, we separate terms with a maximum number of the 
operators a^ and Oq on the right-hand side of (43). Then, from (43), we obtain 

{hz + eq) ganomil^^) = -^{((^a-qflo | a+)) , + {{a+a+a^ \ a+)),} , (44) 

Using the assumption (1) on the C- number representation of the particle creation a^ and 
annihilation ao operators at a zero momentum in the calculation of Green's functions on the 
right-hand side of (44), we find 

{hz + Eg) g^^^^iq, z) = -nou{q) {g^{q, z) + g^^^^{q, z)] . (45) 

It is evident that relations (42) and (45) form a set of algebraic equations in the functions 
g^{q, z) and ganomil^ ^)i from which it immediately follows that 

where the spectrum oj{q) is defined by relation (33). Thus, the poles of the one-particle 
Green's function g^{q, z) when using the assumption (1) about the C- number representation 
of the particle creation a^ and annihilation Oq operators at a zero momentum coincide with 
the poles of the "density-density" Green's function ^^(g, 2;) (see (34)). 

The above consideration is based on the assumption that it is sufficient to consider the 
terms with a maximum number of the operators a^ and ag in the calculation of Green's 
functions at the temperature close to zero. In this case, it is supposed that the interparticle 
interaction is weak. 

Let us now show that the result (46) can be obtained based on the consideration of the 
terms with a maximum number of the operators a^ and Oq without the use of anomalous 
Green's functions ganomi.'i^ ^)- To this end, we write the equation of motion for the function 
((alqaoao | a+))^ (see (40)), 

{hz + Eq) ((a+qaoao | a+))^ = -— ^^u{k){{{a^_^a'^_^^apaoao - a+qa+_j^apa_kao - a^_^aQa'^^^apa_]^] 

Then, on the right-hand side of (47), we separate the terms with a maximum number of the 
operators aQ and ao, 

{hz + Eg) {{atcfiotto I aq))2 = — — {{{a^a^agaoao \ 0+))^ -F a^a+qaoaoao | a'^))^} ,(48) 
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We now use the previously advanced assumption (29) that the operator of the number of 
"condensate" particles A^o = fl([ao is the C-number, rather than the operators a^ and Cq. 
Then, from (48), it immediately follows that 

(hz + Eq) ((a+qfloao I a'^))z = -nou{q) {< Nq > g^{q,z) + {{at^aoao \ a^))^} , (49) 



or 



(( a+ aoao | a+)), = - ^^ ^^"A i TT < ^o > g^{q, z), (50) 

^ ^ {hz}^ + e{q) + nou{q) 

If we substitute relation (50) into (40), we again obtain expression for the one- 
particle Green's function g^{q,z). 

Thus, to obtain the known results for the spectra of excitations in the degenerate Bose 
gas, there is no need for the assumption (1) about the C-number representation of the 
particle creation Gq and annihilation Oq operators at a zero momentum and no need for 
putting into consideration anomalous Green's functions (41). 

Let us now return to relation (34) for the "density-density" Green's function x^il^z), 
from which it follows the necessity of calculating the one-particle distribution function /J 
for "overcondensate" states. Substituting the obtained expression (46) for the one-particle 
Green's function into (38), we find 

■''^ 2hu{q) \ 2T J 2 ^ ^ 

Let us pay attention that the dependence of the one-particle distribution function fj for 
"overcondensate" states on the order of limit transitions T — )■ and g — )■ (see (36)) again 
follows from (51). At the same time, to solve the problem of the braced term fg /{Nq) in 
relation (34), we can reason as follows. The minimum nonzero wave vector g in a specified 
large but finite volume V is proportional to V^^^^. Therefore, in the thermodynamic limit 
V — )■ oo, A^ — )■ oo, N/V = const, the term /J/(A'o) can be considered to be zero provided 
that hq = {Nq)/V 7^ 0, independently of the order of the limit transitions T — ?> and g — )■ 0. 
At first sight, the above results in the limit of the weak interaction completely solve 
the problem of the description of the weakly non-ideal degenerate Bose gas. However, it 
should be taken into account that, as follows from (23), (50), the number of particles in 
"overcondensate" states 

{Nt) = N-{No) = V [d'qfj (52) 
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is nonzero even at zero temperature. Thus, in the hmit T — )■ 0, the number of particles in the 
condensate (A^o) 7^ ^- Hence, the assumption about the C-number behavior of the operator 
^0 = ciof^o, and even more so the operators a^ and qq separately cannot be considered as 
grounded even in the limit of the weak interaction. 

VI. GAP AND THE SELF-CONSISTENT DISTRIBUTION FUNCTION 

In this regard, let us return to the consideration of equation (40) under the assumption 
that we correctly considered the effects of the weak interparticle interaction. Then, taking 
into account the above consideration, we will not assume the C-number behavior of the 
operator A^o = Qo^Oo and the operators a,} and Oq separately. Then, from the viewpoint 
of the perturbation theory on the interparticle interaction, the second braced term on the 
right-hand side of (40) is the higher-order term in comparison with the first term; therefore, 
it can be neglected in the assumption of our interest. Let us pay attention that, from this 
point of view, the result (34) for the "density-density" Green's function remains valid. 

Then, in the limit of a weak interparticle interaction, from Eq. (l40l) we obtain 

The expression for the spectrum of the one-particle excitations Eg is given by 

Eq = eq + nou{q), (54) 

From Eqs. (!38|) . ( 152|) . (153!) . for the hmit of strong degeneracy T <^ Tq, for the distribution 
function of the overcondensate particles, we obtain 

^^ = exp{EjT)-r (^^) 

Therefore, the function /J is finite for g — )■ 0. Moreover, in the limit of strong degeneracy 
T^O 

fj -^ (56) 

for arbitrary values of q, in contrast to the case of the ideal Bose gas. Therefore, the rep- 
resentation ( |22|) for the one-particle distribution function fp is valid, the initial suggestions 
(l2|).( fT8|) are satisfied and the used procedure of extraction of the main terms is self-consistent. 
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As is easily seen from ( H2|) , (H3|) for T — )■ all particles in the considered approximation placed 
in the "condensate" with the momentum p = 0, in contrast to the Bogolybov's theory of 
weakly non-ideal homogeneous Bose gas (the effects of inhomogeneity will be considered 
in further publication). According to the accepted point of view, interaction between par- 
ticles leads to the quantum depletion of the condensate, which means the appearance of 
the particles with p ^ a.t T = 0. However, this statement is based on the assumption 
about coincidence of the excitation spectra for the one-particle Green function and dynamic 
structure factor. Only in this case the condensate density no(T) can be calculated, based on 
the experimental data for S{q, u) and the corresponding collective excitations in the spirit 
of the papers [22-24]. This coincidence takes place in the standard theory of weakly non- 
ideal Bose gas, which is constructed on the assumption about anomalous averages. This 
assumption is extended to liquid He II, where the approximate scheme for calculating uq 
is used. Therefore, the depletion of the Bose condensate is not an experimental fact, but 
is the consequence of applying of the theoretical concept of the identity of one-particle and 
collective spectra, obtained in the anomalous averages approach for weakly non-ideal Bose 
gas. It was shown above that the theories based on C-numher approximations for operators 
in the state with zero momentum are not self- consistent. 

In the developed theory the depletion is absent and, according to the above argumen- 
tation, there is no contradiction between this result and the existing experiments for the 
dynamical structure factor, since the latter cannot be applied to calculate of the particle 
distribution. In general, in a more elaborated theoretical approach the depletion can ap- 
pear, however, even for such approach, the coincidence of the poles of one-particle Green 
function and the dynamical structure factor seems an exceptional occasion. Therefore, the 
Bose particle distribution should be found independently of the experiments on neutron 
scattering. 

For the self- consistent approach under consideration the gap in the spectrum of one- 
particle excitations appears and is given by 

A = ^,^0 = ^0^(0). (57) 

A value of the gap is completely defined by the density of particles in the "condensate" . The 
existence of the gap permits to extend essentially the applicability of the results obtained 
for T — )■ 0. It is obvious, that in many applications the condition T — )■ is equivalent to the 
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condition T <^ A. 

Taking into account (jH]) we can rewrite Eq. fl56|) for the function x^^Hl^ ^) i^ ^^^ form 

On the basis of fl58|) almost all known results for the thermodynamical functions of the 
degenerate Bose gas can be reproduced, as is done in [9-11]. 

Therefore, in contrast with the approach, based on the C-number approximation for the 
operators «(]' and Oq (or applying the C-number approximation for the operator Nq in the 
direct perturbation theory in the presence of the condensate) we find that the spectra of the 
collective and one-particle excitations are different. Both spectra, as is easily seen, satisfy 
the Landau condition for superfluidity. For the one-particle spectrum the Landau condition 
is satisfied for the transitions between the "condensate" and the "overcondensate" state. 
The Landau condition is, naturally, violated for transitions between the "overcondensate" 
states. 

VII. CONCLUSIONS 

Summarizing the above consideration, we can assert that the calculation of the Green's 
functions for the highly degenerate Bose gas on the basis of the C-number approximation 
for the operator A'"o and the straightforward perturbation theory allows the following con- 
clusions. 

(i) The problem of 1/g^ divergence, which arises for the ideal Bose gas, can be solved. 

(ii) The system has two different branches of excitations, i.e., the single-particle and 
collective ones, both satisfying the Landau condition of superfluidity. 

(iii) The single-particle excitation spectrum contains the gap in the region of small wave 
vectors, associated with the existence of the "condensate". 

(iv) The spectrum of collective excitations corresponds to "phonon-roton" excitations 
observed in the experiments on inelastic neutron scattering [23,24,32]. This spectrum is the 
gapless mode which has no connection with the Goldstone theorem. 

(v) There is no need for anomalous averages (quasi-averages) to describe the degenerate 
Bose gas. 
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In this connection, we emphasize that the Goldstone theorem [25] which justifies the gap- 
less branch of excitations is itself the result of the breaking symmetry assumption. Therefore, 
it cannot be applied to the theory considered above, which is constructed without breaking 
symmetry assumption (e.g., without anomalous averages). In our paper, we showed that the 
breaking symmetry is not necessary (in this connection see also [26-29], where the various 
approaches to avoid the breaking symmetry have been developed) and at the same time the 
Bogolyubov result for the collective phonon-roton mode is completely reproduced. In [30,31], 
the analysis of the applicability of the relativistic Goldstone theorem to the non-relativistic 
statistical theory has been done. 

It should also be mentioned that Hugenholtz and Pines, constructing the diagram tech- 
nique for the degenerated Bose gas in [3], reformulated the problem at the beginning, by 
changing, according to Bogolyubov, the operators a^ and ag by C-numbers. On this way, 
they could apply almost automatically the quantum field theory methods to the study of 
the Bose gas with "condensate". In particular, it was shown that, for the C-number repre- 
sentation of the operators a^ and ao, the gap in excitations associated with the one-particle 
Green's function cannot exist. This result is also associated with the anomalous averages 
assumption. 

In the case at hand, it should be noted that this assumption has no any fundamental 
basis, is not used for the ideal Bose gas and, as we showed above, can be avoided for the non- 
ideal Bose gas. Historically, the assumption of breaking symmetry in the Bose gas theory 
played a crucial role and was very useful, but in fact it is not necessary. 

The proposed theory does not contradict the fundamental physical laws; the known ex- 
periments on the collective mode can be qualitatively described since this mode is the same 
as that in the existing theory based on the breaking symmetry assumption. The Landau 
criterion of superfluidity is also fulfilled. 

Therefore, the interparticle interaction in Bose systems leads not only to the drastic 
difference (in comparison with the ideal Bose gas) in the structure of collective excitations, 
which are described by the "density- density" Green's function, but also to the crucial change 
in the distribution function of single-particle excitations and in the single-particle excitation 
spectrum of for "overcodensate" particles. 

Based on the results obtained, the special diagram technique can be developed, similar 
to [19], but with the use of the C-number approximation for the operator Nq. 
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The principal difference between the results of this study and the results of the 
"traditional" C- number approximation for the operators a^ and Oq (as well as for the op- 
erator A'"o) is the existence of the gap in the spectrum of single-particle excitations. The 
above analysis shows that this gap cannot manifest itself in the experiments on inelastic 
neutron scattering in superfluid helium (see, e.g., [32-34]). However, such possibility cannot 
be excluded [16] in the experiments on Raman light scattering. Moreover, in [35-38], where 
such experiments are described, there is a direct indication of the existence of the gap. 

Probably, the gap for single-particle excitations can manifest itself in the experiments on 
the density profile in trapped Bose gases [39]; however, this question cannot be considered 
within the present paper and has to be studied separately. 
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